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Abstract. Quantum ergodicity of classically chaotic systems has been studied extensively 
both theoretically and experimentally, in mathematics, and in physics. Despite this long 
tradition we are able to present a new rigorous result using only elementary calculus. In the 
case of the famous Bunimovich billiard table shown in Fig.l we prove that the wave functions 
have to spread into any neighbourhood of the wings. 



The quantum/classical correspondence is a puzzling issue that has been with us since the 
advent of quantum mechanics hundred years ago. Many aspects of it go back to the New- 
ton/Huyghens debate over the wave vs. corpuscular theories of light. 

On the surface of our existence we live in a world governed by laws of classical physics. That 
does not mean that we know precisely how fluids flow or solids move. They are described by 
highly non-linear rules which are hard to unravel mathematically. Even the simplest classical 
motion, that of a ball bouncing elastically from confining walls poses many unanswered questions 
- see http: //www. dynamical systems.org/billiard/ for a fun introduction. 

If we investigate deeper, or if we simply use any modern technical device, we come in contact 
with quantum mechanics. It is governed by a different set of rules which mix wave and matter. 
The simplest description of a wave comes from solving the Helmholtz equation: 

(-A-A 2 )u = 0, A = d^ + dy, (as,y)efi, u\ dn =0. 

Here we put our wave inside of a two dimensional region Q. In classical wave mechanics the 
limit A — > oo is described using geometrical optics where the waves propagate along straight 
lines reflecting in the boundary dfl. Roughly speaking, we expect something similar in the clas- 
sical/quantum correspondence with the Helmholtz equation replaced by its quantum mechanical 
version, the Schrodinger equation. For many fascinating illustrations of this we refer to the web 
art gallery of Rick Heller: http://www.ericjhellergallery.com 

Many researchers on different aspects of semiclassical analysis have been interested in the the 
correspondence of solutions to the equation above and the classical geometry of balls bouncing 
from the walls of f2: Backer, Cvitanovic, Eckhardt, Gaspard, Heller, Sridhar, in physics, and 
Colin de Verdiere, Melrose, Sjostrand, Zelditch, in mathematics, to mention some (see [2 ,[3J for 
references to the physics literature, and 0,0, |6| for mathematics). 

Billiard tables for which the motion is chaotic are a particularly interesting model to study 1, . 
One of the most famous is the Bunimovich billiard table shown in Fig.l. By adding two circular 
"wings" to a rectangular table the motion of a reflecting billiard ball becomes chaotic, or more 
precisely, hyperbolic, in the sense that changes in initial conditions lead to exponentially large 
changes in motion as time goes on. 



tOf course one would not want to play billiards on a table like that, and a completely integrable rectangular 
one can pose enough of a challenge. While discussing billiards and the classical/quantum correspondence we 
cannot resist mentioning that Pyotr Kapitsa (Nobel Prize in Physics '78) was fond of saying that trying to detect 
the quantum nature of physical processes at room temperature was like trying to investigate the physical laws 
governing the collision of billiard balls on a table aboard a ship going through rough seas. 
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Figure 1. An experimental image of bouncing ball modes in a Bunimovich 
stadium cavity - see pQ and http://www.bath.ac.uk/~pyscmd/acoustics. With 
a certain amount imagination one can see our theorem in this picture. 

As a model for studying quantum phenomena in chaotic systems this billiard table has become 
popular in experimental physics. A genuinely quantum example is shown in Fig. 2 - it comes from 
the scanning tunnelling microscope work of Eigler, Crommie, and others 




Figure 2. Quantum corral in the shape of the Bunimovich stadium. Courtesy 
of IBM® Research. 

One question which is still mysterious to mathematicians and physicists alike is if the states 
of this system (that is, solutions of the equation above) can concentrate on the highly ustable 
closed orbits of the classical billiard. Quantum unique ergodicity states that there is no such 
concentration - see [8 ,[9|, 6 and references given there. In the arithmetic case, that is for 
billiards given by arithmetic surfaces where the motion is given by the geodesic flow, spectacular 
advances have been recently achieved by Bourgain, Lindenstrauss ^H], and Sarnak, while for the 
popular quantization of the Arnold cat map impressive results were produced by Bonechi, De 
Bievre, Faure, and Nonnenmacher and also by Kurlberg and Rudnick |12j . 

Here we describe an elementary but striking result in the billiard case. It follows from adapting 
the first author's earlier work in control theory. Although motivated by the more general aspects 
of [13 we give a simple self contained proof. 
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For a state u(x,y) what counts is its probability density, \u{x,y)\ 2 dxdy - we assume here that 
\u\ 2 is normalized to have integral 1 over f2. We say that it is bounded from below in a region if 
its integral over that region is bounded from below by a positive constant. With this terminology 
we have, roughly speaking, 

Theorem. For any normalized state of the Bunimovich billiard table, the probability density in 
any neighbourhood of the wings of the table is bounded from below independently of the energy 
A 2 . 

In particular, the result says that single bouncing ball orbits (that is orbits following an interval 
perpendicular to the horizontal straight boundaries) cannot produce localized waves. Our result 
allows concentration on the full invariant set of all vertical orbits over R - that is consistent 
with the existing physical literature - both numerical and experimental - see f2] and pp. In ^3] 
we show a stronger result, namely that the neighbourhood of the wings can be replaced by any 
neighbourhood of the vertical intervals between the wings and the rectangular part. The proof 
of that predictable (to experts) improvement is however no longer elementary and is based on 

The proof of theorem depends on the following unpublished result of the first author (see ^2] 
for detailed references and background material): 

Proposition. Let A = d 2 + d 2 , be the Laplace operator on the rectangle R — [0, l] x x [0, a] y . 
Then for any open u> C R of the form lo = uj x x [0, a] y , there exists C such that for any solutions 
of 

(-A - X 2 )u = f + d x g onR, u\ dR =0, 
with an arbitrary X > we have 

[ \u(x,y)\ 2 dxdy<c([ (\f(x,y) 2 \ + \g(x,y)\ 2 )dxdy + [ \u(x,y)\ 2 dxdy) . 
Jr. \Jr Juj J 



Proof. We decompose u, and / + d x g in terms of the basis of L 2 ([0 7 a}) formed by the Dirichlet 
eigenfunctions ek(y) = y/2/a sm(2kiry / a) , 

(1) u(x,y) = ^2e k (y)u k (x), f(x, y) + d x g{x, y) = ^ e k (y)(f k {x) + d x g k (x)) 

k k 

we get for u k , f k the equation 

+ z)u k = fk + d x 9k, Ufc(0) = Ufc(l) = 0, z = X 2 - (2kn/a) 2 . 
It is now easy to see that 

(2) / \u k (x)\ 2 dx<c( f (\f k (x)\ 2 + \g k (x)\ 2 )dx+ f |ti fc (x)| 2 da:V 

JO \Jo Jul* / 

where C is independent of Ai. In fact, let us first assume that uo x — (0,(5), 8 > 0, and z = A 2 , 
with ImAi < C. We then choose \ € C£°([0, 1]) identically zero near and identically one on 
[<5/2,l]. Then 

(pi + A 2 ) (x«fc) = F k , F k = x{fk + d x9k ) + 2d xX d x u k + d 2 X u k . 
We can now use the explicit solution given by 

1 r 

X(x)u k (x) = — / sin(Ai(a; - y))F k (y)dy . 



N. BURQ AND M. ZWORSKI 



All the terms with d x g k and d x u k can be converted to gk and u k by integration by parts (with 
boundary terms at both ends). Due to the X^ 1 factor that produces no loss and the estimate 
follows. The argument is symmetric under the x i— > —x change, so we can place our control 
interval anywhere. 

It remains to discuss the case z < —C < 0. Then the estimate (2J follows from integration by 
parts (where now we do not need u> x ): 



fk(x)u k (x) - g k {x)d x u k (x) dx = (fk(x) + d x g k {x))u k (x)dx 



(-d 2 x - z)u k (x)u k (x)dx — I (\d x u k (x)\ 2 + \z\\u k (x)\ 2 ) dx . 

Jo 

By the Cauchy-Schwartz inequality, the left hand side is bounded from above by 

1 ~ l 

{\fk(x)\ 2 + \g k (x)\ 2 ) dx^j 2 Qf (\u k (x)\ 2 + \d x u k (x)\ 2 ) dx 

Since \z\ > C > 0, (2J follows from elementary inequalities (see Lemma 4.1] for a general 
microlocal argument). We can now sum the estimate in k to obtain the proposition. □ 

We can now present a more precise version of the theorem. For a yet finer version we refer 
the reader to [T^l Theorem 3'] and [H3 Fig.5]. 

Theorem'. Consider fl the Bunimovich stadium associated to a rectangle R. With the convention 
of Fig. 1, let R± be any rectangle with the horizontal sides contained in the sides of R, strictly 
contained in R, and with R\R± having two components. 

There exists a constant C depending only on fl and R\ such that for any solution of the 
equation 

(-A-\ 2 )v = f, u\ dn =0, A>0, 

we have 

\v{x,y)\ 2 dxdy < C\ \ \f(x,y)\ 2 dxdy+ / \v{x,y)\ 2 dxdy) . 

\Jn Jn\R! ) 

The "wings" of the billiard table in the original statement are given by O \ R% , We apply the 
second theorem with / = to obtain the first one. 

Proof. Let us take x, y as the coordinates on the stadium, so that x is the horizontal direction, 
y the vertical direction, and the internal rectangle is [0, l] x x [0,a] y . Let us then consider u and 
/ satisfying (—A — \ 2 )u = /, u = on the boundary of the stadium, and x( x ) £ C£°(0, 1) equal 
to 1 on [e, 1 — e]. Then x(x)u(x,y) is solution of 

(-A-A 2 ) X u = x/+[A,x]" ini? 

with Dirichlet boundary conditions on OR. Since [A,%]u = 2d x (x' u ) ~ x" u we can a PPly the 
proposition to obtain 

\x(x)u(x,y)\ 2 dxdy < C ( / \x(x)f(x, y)\ 2 dxdy + / \u(x , y)\ 2 dxdy 
\Jr Juj e 

where u E is a neighbourhood of the support of d x x- Since we can choose it to be contained in 
R\Ri, the theorem follows. □ 
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We conclude by remarking that the same argument holds in the setting discussed recently in 
|15j and 0, since in the argument above the rectangle can be replaced by a torus. 
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